import numpy as np
import sympy as sy
al=np.array([[1,2,3,-4]]).T
a2=np.array([[0,1,-1,-1]]).T
a3=np.array([[1,3,0,1]]).T
a4=np.array([[0,-7,3,1]]).T
#按水平方向(列顺序)堆叠数组构成一个新的数组
A=np.hstack((al,a2,a3,a4))
C=sy.Matrix(A)
print('矩阵A为:\n',A)
print('矩阵A的秩为: ',np.linalg.matrix_rank(A))
print('矩阵A的行列式为:',np.linalg.det(A))
print('矩阵A的行最简形为:\n',C.rref())
